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Abstract 



We provide evidence for the influence of thermal fluctuations of Goldstone modes on the chiral 
condensate at finite temperature. We show that at fixed temperature, T < T c , in the vicinity of 
the chiral transition temperature this leads to a characteristic dependence of the chiral conden- 
sate on the square root of the light quark mass (mi), which is expected for 3-dimensional models 
with broken O(N) symmetry. As a consequence the chiral susceptibility shows a strong quark 
mass dependence for all temperatures below T c and diverges like l/^/m7 in the chiral limit. 



Key words: QCD thermodynamics, chiral symmetry breaking, Goldstone modes 
PACS: 11.15.Ha, ll.10.Wx, 12.38.Gc, 12.38.Mh 



1. Introduction 

Establishing the properties of the chiral phase transition in QCD is one of the outstand- 
ing problems in lattice simulations of QCD thermodynamics. Although several attempts 
have been undertaken to establish the unique scaling behavior in the limit of vanishing 
light quark masses no clear-cut evidence for the expected O(N) scaling, related to the 
underlying chiral symmetry of the QCD Lagrangian, has been found. This may have 
several reasons related to common problems in lattice calculations: cut-off effects, finite 
volume effects, explicit flavor symmetry breaking and/or too large quark masses. 

We will discuss here yet another effect that influences the scaling behavior of ther- 
modynamic quantities in the vicinity of the finite temperature chiral transition of QCD 
and may compete with the universal scaling in the vicinity of T c . This is related to sin- 
gularities in chiral observables, e.g. the chiral susceptibility, that are induced at ANY 
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Fig. 1. Chiral condensate calculated on lattices of size ATj X 4 in (2 + l)-flavor QCD for various values 
of the light quark mass (mi) in units of the strange quark mass (m s ). N& from 32 for the two lightest 
quark mass values to 8 for the heaviest masses. For small quark masses the critical coupling for the chiral 
transition temperature on this size lattices has been determined to be in the range j3 c = 3.305 — 3.31 [7]. 

temperature in the chiral symmetry broken phase. The existence of these singularities 
induced by fluctuations of Goldstone modes are well known from the analysis of 3 and 4 
dimensional statistical models and field theories with global O(N) symmetries that are 
spontaneously broken at low temperature [1—4]. 

2. Chiral condensate and susceptibility for T<T C 

At vanishing temperature properties of QCD in the light quark sector are described 
by an effective 4-d, 0(4) symmetric Lagrangian; chiral symmetry is broken and, what is 
of interest here, at non-zero quark mass the light quark chiral condensate, {iptf)), receives 
logarithmic corrections proportional to to/ In to; that are due to the presence of light 
Goldstone modes [2]. This leads to a logarithmic divergence of the chiral susceptibility, 
Xm ~ d(ipij}} /dm. In three dimensions the analog to these logarithmic terms are correc- 
tions that are proportional to the square root of the quark mass [1]. As a consequence 
Xm will diverge in the chiral limit like Xm ~ 1/^/to,;. 

In the vicinity of the QCD phase transition a 3-d, 0(4) symmetric effective theory is 
expected to describe the universal properties of the transition [5] which will also control 
the properties of the chiral condensate for temperatures close to T c . We thus expect that 
in the chiral limit Xm will diverge for all values of the temperature T < T c , 



where 1/8 — 1 ~ 0.79 for 3-d 0(4) as well as 0(2) symmetric spin models. 

The appearance of a square root singularity has, in fact, been established in the vicinity 
of the chiral phase transition of an SU(3) gauge theory coupled to adjoint fermions [6] 
where it is apparent already for moderately small values of the quark mass. Similar 
evidence for the role of light Goldstone modes in QCD, however, has been missing so far. 

In Fig. 1 we show the light quark chiral condensate, {'tpip) calculated in (2+l)-flavor 
QCD with improved staggered fermions (p4-action) on lattices of size N% x 4 for various 
values of the gauge coupling. The figures show results for (ipi/)}, plotted versus the light 
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Fig. 2. Disconnected part of the chiral susceptibility calculated on lattices of size X 4 in (2 + l)-flavor 
QCD for various values of the light quark mass. 

quark mass (mi) expressed in units of the strange quark mass (m s ) (left) and the square 
root of this ratio (right). The right hand figure indeed suggests that the strong curvature 
of (iptp) as function of mi is well accounted for by a yjmi dependence for T < T c . In 
any case, the weaker than linear dependence of (tpip) on mi leads to a strong quark mass 
dependence of the chiral susceptibility, Xm- This susceptibility receives two contributions, 
usually referred to as disconnected and connected parts, X m = Xdis + Xcon, with 

((TrZY 1 ) 2 ) - (Tr^ 1 )) , ^ = H^R'f"^)) ■ 

Here Di denotes the staggered fermion matrix for light quarks. In general, both contri- 
butions are sensitive to the Goldstone modes and will receive divergent contributions, 
Xdis.con ~ 1/y/mi [4]. We will restrict our analysis here to the disconnected part. 

In Fig. 2 we show Xdis/T 2 calculated for various values of mi. In these calculations the 
strange quark mass is close to its physical value and mi reaches values well below the 
physical one; a ratio mi/m s ~ 0.04 corresponds in the continuum limit to the physical 
quark mass ratio. The strong quark mass dependence in a wide temperature interval 
is obvious from Fig. 2 (left); its dominant contribution in the low temperature regime 
is indeed well described by x-m ~ ^l\/mi as can be seen from the rescaled data shown 
in Fig. 2 (right). We note that this seems to hold for all quark masses, mi < 0.1m s , 
although at present our data sample for the lightest quark mass is not yet good enough 
to substantiate this statement also for this light mass value, which corresponds to a light 
pseudo-scalar mass of about 70 MeV. We furthermore, note that the square root behavior 
also dominates the quark mass dependence in the peak region, i.e. there is no evidence for 
a stronger singularity, which is expected to arise from the singular part of the partition 
function of 3-d, 0(4) (or 0(2)) symmetric models, x s ™av,la,r _ TO -°- 79 . 



Xdis _ N T / 
T 2 ~ l&Nl I 



3. Towards the continuum limit of (2+l)-flavor QCD 



The general features of the quark mass dependence of the disconnected part of the 
chiral susceptibility are also apparent in calculations currently performed by the hotQCD 
collaboration on lattice with temporal extent N T = 8 [8], i.e. closer to the continuum 
limit. In Fig. 3 we show Xdis/T 2 and the rescaled susceptibility. The left hand part of the 
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Fig. 3. Disconnected part of the chiral susceptibility calculated on lattices of size 32 3 x8 in (2 + l)-flavor 
QCD for mi/m s = 0.05, 0.1 and 0.2, respectively. 

figure shows that for temperatures above the transition region the susceptibilities quickly 
become quark mass independent. This is expected as the chiral condensates will be linear 
in the quark mass in the chirally symmetric phase of QCD. Below the transition region 
we again observe a strong quark mass dependence that is consistent with ~ l/^/mi. 

The weak quark mass dependence of Xm in the high temperature regime suggest that 
the susceptibility can already be described by the limiting form for the mi — > limit, 

oo , T <T C 

A consistent picture for the chiral susceptibility thus seems to emerge despite the missing 
evidence for O(N) scaling for its peak height. 
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